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Abstract 

The search for an understanding of an energy source great enough 
to explain the gamma-ray burst (GRB) phenomena has attracted much 
attention from the astrophysical community since its discovery. In this 
paper we extend the work of K. Asano and T. Fukuyama, and J. D. 
Salmonson and J. R. Wilson, and analyze the off-axis contributions to the 
energy-momentum deposition rate (MDR) from the vv collisions above 
a rotating black hole/thin accretion disk system. Our calculations are 
performed by imaging the accretion disk at a specified observer using 
the full geodesic equations, and calculating the cumulative MDR from 
the scattering of all pairs of neutrinos and anti-neutrinos arriving at the 
observer. Our results shed light on the beaming efficiency of GRB models 
of this kind. Although we confirm Asano and Fukuyama's conjecture as 
to the constancy of the beaming for small angles away from the axis; 
nevertheless, we find the dominant contribution to the MDR comes from 
near the surface of the disk with a tilt of approximately 7r/4 in the direction 
of the disk's rotation. We find that the MDR at large radii is directed 
outward in a conic section centered around the symmetry axis and is 
larger, by a factor of 10 to 20, than the on-axis values. By including 
this off-axis disk source, we find a linear dependence of the MDR on the 
black hole angular momentum (a). In addition, we find that scattering 
is directed back onto the black hole in regions just above the horizon of 
the black hole. This gravitational "in scatter" may provide an observable 
high energy signature of the central engine, or at least another channel 
for accretion. 



'Department of Mathematics, North Carolina State University, Raleigh, NC 27695-8205 
t Computer and Computational Sciences Division, Los Alamos National Laboratory, Los 
Alamos, NM 87545 

* Theoretical Astrophysics Group (T-6, MS B288), Theoretical Division, Los Alamos Na- 
tional Laboratory, Los Alamos, NM, 87544, and Air Force Weapons Laboratory, High Power 
Microwave Division, Kirtland AFB, NM 87117; Manuscript Contact E-mail: wam@lanl.gov 



1 



1 Neutrino Scattering About a Low Mass Ac- 
creting Kerr Black Hole 

The search for an understanding of an energy source great enough to explain 
the gamma-ray burst (GRB) phenomena has attracted much attention from the 
astrophysical community since its discovery. One of the more recent candidates 
for the GRB model is based on the fire-ball model, [jj], |], ^] The proposal for 
the central engine for such bursts concerns vv production and their subsequent 
scattering in the environs of a dense accretion disk/low-mass black hole system 
[^[ [| . Such a system may be formed by the merger of a black hole and a neutron 
star, the merger of two neutron stars, by a collapsar, or by a supernova |^|, [?], ||, ^| . 
The astrophysical details of the geometry or environment of such systems are 
currently hidden from us both observationally and computationally, although 
this situation may change in the near future with the development of gravity 
wave observatories and more sophisticated general relativistic astrophysical sim- 
ulation capabilities. For the mean time in this paper we focus our attention on 
a rather simplified model for the scattering of neutrinos and anti-neutrinos in 
the vicinity of a low-mass rotating black hole. Our goal is to extend the recent 
work of Asano & Fukuyama (AF) (h], |ll| and Salmonson and Wilson Q (SW) 
by providing the first general-relativistically covariant off-axis calculation of the 
energy momentum deposition rate (MDR). We focus our attention, as did the 
previous authors, on the general relativistic enhancement or degradation of the 
deposited MDR. 

Since we now intend to calculate the off-axis contributions, the direction 
of the energy-momentum 4-vector can no longer be assumed, as opposed to 
the along axis case, where symmetry forces the net energy upward along the 
symmetry axis. An accurate treatment of off-axis contributions requires a full 4- 
dimensional spacetime calculation, as opposed to the analytic and semi-analytic 
calculations performed by AF and SW. Therefore, we derive the full energy- 
momentum (hereafter referred to as " momenergy" ) [[l3| deposited in the vicinity 
of the black hole/accretion disk system. The calculations are based on the basic 
principles as discussed in Misner, Thorne and Wheeler. |l4jl We consider the 
direction of the momenergy in the frame of an observer at a given point above 
the disk. 

Our aim in this work is not solely to study and draw conclusions on the 
relativistic case of neutrino and anti-neutrino collisions above the accretion disk 
for GRB models, but also to develop a method to provide a more in-depth 
understanding of the scattering mechanisms in the environs of small accreting 
black hole systems. For the above purpose, we consider an idealized situation 
which includes gravitational redshift, the bending of neutrino trajectories, and 
the redshift due to accretion disk rotation. 

This paper is organized as follows: In Sec. ^ and Sec. || we define the ap- 
propriate global spacetime metric and establish the frame of the observer. The 
derivation of the momenergy deposition rate (MDR) is discussed in Sec. f|. Sec. |^ 
deals with the computation of the MDR in an observer's frame with regards to 
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evaluating the number densities and the generation of the 4-momentum from 
the v-v scattering. In Sec. ^| we give an estimate of the cumulative effect of the 
scattering on the energy deposition at a radius of 50M from the black hole. Fi- 
nally, in Sec. ^we draw some conclusions, discuss the applicability of our results 
and methods on a larger scale, and make recommendations for future work in 
this area. 



2 The Global Spacetime Metric 

We use the Boyer-Lindquist (BL) spacetime metric in our neutrino ray tracing 
code. 

ds 2 = g tt dt 2 + 2 g t( /, dt dip + g rr dr 2 + g gs dO 2 + g^ dp 2 (1) 

where 

2Mr 

9tt = + ( 2 ) 

2Mr 2 . . 

9t4 = —asm 9 (3) 

E 

9rr — ^ 

gee = E (5) 

94>4> = a 2 sin 2 9 + r 2 + a 2 ^j sin 2 6 (6) 

E = r 2 + a 2 cos(9 (7) 

A = r 2 + a 2 - 2Mr (8) 



(4) 



with < a < M. 



3 The Frame of the Observer 

The momenergy deposition rate (MDR) is computed in the frame of an observer 
resting with respect to the global coordinates it, r, 9, 4>) at a point with coordi- 
nates r,9,(f> = const. Spacetime is stationary which implies that spatial pictures 
are sufficient. Also, <f> can be fixed (say cj> = 0); however, the observational ef- 
fects are nonetheless ordinarily only representable in 3-dimensions because of 
frame dragging. 

The frame of the observer is defined by four orthonormal basis vectors 
({e^}? =0 ). We define these basis vectors such that (1) the basis is a Lorentz 
frame (e^ • = ??^), and (2) the 4-velocity, or temporal basis vector of the 
observer (u = eg) is parallel to the time component of the BL coordinate frame, 
e = d t . 

In the global coordinate frame 

e 6 = ad t (9) 
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eg • e 6 = a 2 d t ■ d t = .goo" 2 = -1 ==> a = (10) 



This leads to 



Likewise 



' ft (11) 
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Vector is determined by the conditions = a 9*4-/3 d^, e^-dt = 0, e^-e^ = 1, 
which leads to 

*i = J m 2 (-^d t +d,) (14) 

4 The Momenergy Deposition Rate in the Frame 
of the Observer 

The MDR in the frame of the observer is given by 

d?p v d 3 p p 



MDR =11 fu{Pu, ■ ■ ■) U(Pp, ■ ■ ■) i>|v„ - v v \s v e v } {p v + p p ) 



(15) 

where /„, fp are number densities in phase space (Lorentz invariant), and a 
is the rest frame cross section. The expression in curly brackets is Lorentz 
invariant and can be computed as 

{cr|v„ - V b \e„£ d } = (-£ v £y + Pv ■ Pvf (16) 

where 

G% = 5.29 x lCT 44 cm 2 MeV~ 2 (17) 

D = l±4sin 2 6>vi/ + 8sm i 9 w (18) 

In Eq. (J18|), 0w is the Weinberg angle, sin 2 9w = 0.23 and the plus sign is used 
for v e v e pairs, while minus is used for and v r v T pairs. 

In order to perform integration in spherical coordinates, we use expressions 
Pv = £v and d p u = e v de v dQ, v , where £l v is the unit vector in the direction 
of p y and <ifi„ is the solid angle element. Equation (|l5|) then reduces to 

DG 2 f f 

MDR = — — ^ / J v f u {0. v -^l iJ -l) 2 {e l ,n i , + e i ,0.p)elelde v de v dn v dn v 



3tt 



(19) 
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5 Computing the Number Densities f„, f D 




Figure 1: The geometry behind the computation of the total momenergy depo- 
sition rate in the observer's orthonormal frame (hatted basis vectors). 

The number densities /„, fp are Lorentz invariant and conserved along the 
v, v world lines (null geodesies). If p = (e, p) is the neutrino (anti-neutrino) 
4-momentum in the observer frame (e = |p|), and = (£i|,pi|) is the 4- 
momcntum of the same neutrino in the comoving frame of the disk at the point 
of emission {e\\ = — u e • pp) then /„ is computed as 

s.-i^-)-^{j f )-U' fb '"* 1 T (20 > 

where l v is the specific intensity of radiation at a given frequency v. Note that 
the ratio ij- is Lorentz invariant by Liouville's Theorem. The value of £||„ in this 
equation is determined by the 4 momentum of the neutrino parallel translated 
to the disk along the world line of the neutrino p^ and on the 4-velocity of the 
particle in the disk emitting the neutrino u e . The latter is computed as 
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where pi 



d A = n = Ml/2 m) 

dt r^ + aM 1 / 2 v ' 



for the direct circular orbit of the emitting particle. 
The normalization condition 

Urf-u d = -l (23) 

written as 

(dt \ ^ 
— J (g tt + 2g^n + g^n 2 ) = -1 (24) 

allows the evaluation of it 4 as 

«* = ^ = [- Ut + 2 5t0 n + 500 « 2 )] ~* (25) 
The value of Ud is given by 

u i =/« t J 0,0,u*^\ (26) 



The steps for evaluating /j,(p) are as follows: 

1. p — ► p„ = in the frame of observer (supplied by the integrating 
program) 

2. p„ = (p* , p r , p e , p^) Transformation to global coordinates 

3. Tracing geodesies starting at (t, r, 8, 4>) bs tangent to p„ to the disk. Pick 
the geodesic parameter A such that -4r = p„, integrate back to 9 = 5 (A = 
at the observer, A < at the disk). This procedure yields (t, r, 8, <f)disk 
and p||„ = {d/d\) dls k = (p* , p r , p 9 , p^) disk 

4. Computing u(t, r, 0, = u d 

5. Computing em = — • p y 



6. Computing f(p) = f(e\\) using equation ( 20 ) 

Assuming the emission rates of the neutrinos and the anti-neutrinos at the 
disk are identical, and assuming the neutrinos are emitted isotropically as a 
Planck black-body spectrum characterized by the temperature Tdisk hi the frame 
comoving with the disk material, then the spectrum of neutrinos that our ob- 
server sees is also a Planck spectrum. However, the temperature characterizing 
our observer's spectrum is scaled by the frequency shift of the neutrinos. 

To = —Tdisk (27) 
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Following the derivation of Salmonson and Wilson) 12 , the MDR integral (Eq.|l9|) 
can be transformed into a 4-dimensional integral by analytically integrating over 
£„ and e p . 

MDR = 77t3 (^ G f J J (kTo)\kTof(£l v ■ - lftl v dQ v (28) 

Here, the temperature Tq is the temperature of the black body spectrum of the 
neutrinos as measured by the observer O. It will ordinarily depend on both the 
sky angles (ip, £) we integrate over. In a similar way Tq = Ta{tj},^). Using 
Eqn. [27j we can rewrite the momenergy deposition rate in terms of the disk 
temperature. 

MDR = aI(r o ,0 o ) (29) 

where tq and do are the BL coordinate location of the observer, the integral 
function is driven by the redshift of the neutrinos and anti-neutrinos arriving at 
the observer, 

I{ro,0 o ) = J J (jf\ (*-f) {n v -a s -\fa v dSl v dSl B (30) 

V V 

and the constant a depends strongly on the temperature of the disk as measured 
by the Keplerian observer. 

7ir 3 ((5)DG 2 F 9 
" = ( kT ^sk) 

= 1.84 xlO 32 ( kT f k r X er9S , (31) 
\10MeV ) sec- cm? y ' 

In order to calculate the generation of 4-momcntum from the v-T> scattering, 
we must first image the accretion disk in the observer's reference frame. We do 
this by using a full 3-dimensional ray-tracing code for a Kerr spacetime. The 
input needed by the ray-tracing code in order to image a given pixel of the 
disk is the location of the observer [O = {t, r, 9, 4>}) and the components of the 
4-momentum (p) in Boyer-Lindquist coordinates. 



P = P°e„ + P r e? + p 8 e § + p^e^ 
= p°e + p r e r + p e e e + p% 

We can express the 4-momentum of the neutrino in the observer's frame by 
using the two spherical sky angles (ip and £) shown in Fig. p|. 

P 6 = \P\ - e v 
p r = \p\ cos (tp) 

p = \p\ sin (if)) cos (^) 

p^ = \p\ sin (ip) sin (£) 
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Figure 2: Spherical coordinate sky angles (ip and £) used, in part, to generate 
the 4-momentum p in the observers rest frame. 



The magnitude of the 3-momentum and the time component of the 4- 

momcntum, were obtained by assigning an energy to the neutrino in the ob- 
servers frame (e v ) and using the normalization condition, p • p = 0. We 
can recover these components in the Boyer-Lindquist coordinates by inverting 
Eqnfl.(0@g|0). 



fjti 



P 
P° 



V~9tt y g^gtt - 9i 4 \ gtt 
cos (ip) 



— I e v sin (ip) sin (£_ ) 



\J 9rr 

e 



— sin (ifj) cos (^) 
gee 



— — £ v sin W sin (0 

5005tt - g t< f, 



To reconstruct the image of the entire disk at O, we integrate over all values 
of the sky angles (ip € {0, 7r} and £ € {0, 27r}). We use the ray-tracing codeJi6[ 
to calculate the 4-momentum of the neutrino (p||)when it hits the disk. 

In Fig. H we provide an illustrative example of the imaging of an accretion 
disk extending from r w 4.233M (the inner-most stable orbit of a Kerr black 
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hole of mass M and specific angular momentum a = 0.5) out to r = 10M. 
Here the observer O is located at BL coordinates r = 8M and = 3ir/8. The 
image of the disk is color coded as to the energy of the neutrinos reaching the 
observer. This image is used to calculate the MDR (Eqn. ^). We also provide 
a convergence plot for the MDR in Fig. ||. 

6 Estimate of the Transport of the Momenergy 
out to the Jet-Producing Region 

Up to this point we have calculated the total momenergy density created in 
a given observer's orthonormal frame per unit 3-volume per unit proper time. 
These results are posed properly in general relativity as they represent a mean- 
ingful local observation. However, for the GRB application we address in this 
paper it would be useful to determine the electron distribution function for an 
observer at large radii (e.g. r = 50M) from the accretion disk. This location 
would approach the inner envelope of a hypothetical stellar shell surrounding 
the black- hole/accretion disk system. At such a location we could construct a 
source term for the fluid at the inner envelope of the star that would provide 
the energy for a jet which would ultimately punch through the outer atmo- 
sphere of the star - a jet that would be the generator of a fireball-type solution. 
Unfortunately, the construction of such a source is a computationally difficult 
problem even under the ideal assumptions that the electrons created evolve as 
a collisionless fluid and are highly energetic (E 3> m e ). If we want to know how 
many electrons with momentum p arrive at an observer located at spacetime 
event x, then we would need to construct the geodesic into the past tangent to 
p at x and sum up the electron creation functions along and tangent to that 
geodesic. fDJ This would require modifying the calculation of the MDR (Eq. |l5|) 
to calculate the probability distribution functions for the creation of electrons at 
each point. We would then need to construct and store these electron creation 
functions on a dense grid in the z-x plane perpendicular to the disk and inte- 
grate these creation distribution functions over the 3-parameter (2 sky angles 
and electron energy defining p) family of geodesies of the observer at x. Al- 
ternatively, we could solve a 3-dimensional collisionless Boltzmann equation for 
the Kerr spacetime geometry. [|l8) Unfortunately, both of these approaches, while 
correct in a GR sense, are computationally demanding in the extreme. While 
there is no technical reason that would prevent us from carrying out these cal- 
culations, practical considerations prevent us from completing them in the near 
future. Consequently, we elect to pursue a simpler computational strategy for 
the estimation of the GRB jet source term. 

The difficulties with a general covariant treatment originate with the Kerr 
black hole and not with the Schwarzschild solutions. In particular, we currently 
calculate the total 4-momentum deposited per unit 4-volume in the observer's 
orthonormal frame. Any attempt to assign a 4-volume to the observer in BL 
coordinates (or any coordinate) would break general covariance. In other words, 
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Figure 3: A color rendering of the sky map of an observer located above a ro- 
tating (a — 0.5) black hole of mass (M). The observer is located a distance 8M 
from black hole at an inclination 9 = 3n/8. The thin Keplerian accretion disk 
extends from the inner-most stable orbit (r w 4.233M) out to r — 10M. We 
assume that the neutrinos are emitted from an isothermal disk of temperature 
Tdisk = 1 (as measured by the local comoving Keplerian observers). The spec- 
trum of the neutrinos and anti-neutrinos are assumed to be a Planck spectrum 
characterized by this temperature. The long thin yellow-colored coordinate axis 
represents the observer's orthonormal f coordinate axis. The thicker red arrow 
originating from the center of the observer's celestial sphere represents the sum 
of all the scattered momenergy as projected on the observer's spatial frame. The 
image of the accretion disk is shown color coded by the energy of the arriving 
neutrinos as measured in the observer's frame. In particular, the thin white, 
saturated band on the disk image represents zero frequency shift. The neutrinos 
from the blue region (darker region on the back of the sphere) are blue shifted 
due to the closing relative velocity of the observer and the rotating disk mate- 
rial, and the reddened regions indicate the gravitational red shift effect from the 
black hole. The non-axial symmetry of the disk image, black hole and region 
between the horizon and the inner edge of the disk is due to the black hole's 
rotation (i.e. Lense-Thirring effect, frame dragging or gravitomagnetism). This 
image was generated using a grid on the sphere of dimensions 800 x 1600. 
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Figure 4: A demonstration of the convergence of the energy deposition rate 
and tilt of the scattered electrons/positrons as a function of the number of pixels 
used to image the accretion disk. In this particular convergence test, we use 
an observer located at r = 8M and 9 = 7vt/16 in an (a = 0) Schwarzschild 
spacetime. We chose the location of this observer to correspond to the region 
of maximal momenergy production (as shown in the next figure). We plot 
the scattering integral / in the left plot, and the angle ip in the right plot 
as a function of the number of grid points used to resolve the accretion disk. 
The maximum resolution used in this convergence test involved tracing 2 x 10 6 
geodesies over the entire celestial sphere. Of these geodesies, approximately one 
in seven, or 3 x 10 5 , hit the disk. 
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Figure 5: Pictorial representation of our simulations for three values of the 
spin of the black hole. The first row represents a Schwarzschild black hole 
(a = 0), while the third row an extreme Kerr hole (a = 1). The second row 
represents our results for a moderately rotating black hole (a = 0.5). The 
accretion disk in each of these cases extend from 10M in toward the black hole 
down to the inner- most stable orbit. The first column shows a contour plot of 
the magnitude of the time component of the momenergy deposition rate in Kerr 
coordinates, the second column depicts a contour plot of the angle, or tilt, the 
spatial components (in Boyer-Lindquist coordinates) of the MDR make with 
the symmetry z-axis. In the third column, we combine the information in the 
first two columns and pictorially show the spatial Boyer-Lindquist components 
of the momenergy deposition rate for various locations about the black hole 
accretion disk system. One can observe the dominant contribution to the MDR 
above the surface of the disk produced by the vv collisions as well as their pitch 
angle. Values of the MDR can be obtained by multiplying the numbers in the 
graph by a and by the cube of the black hole's mass. 
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the spatial slices of the BL metric are not orthogonal to the time lines (g^t ^ 0), 
and thus the 3-space of our observers is not integrable. Nevertheless, if we 
take the total MDR calculated in the last section and transform this to the BL 
coordinates and multiply by the "effective BL spatial 3- volume" of our observer, 
we can obtain an approximation to a 4-momentum per proper time of electrons 
created at that location. This is a coordinate-dependent quantity which should 
be adequate for a qualitative understanding as long as the observers are not 
too close the the horizon of a rotating black hole. Moreover, electron-positron 
production close to the horizon should not be expected to contribute significantly 
to the results at large radii, as the electron-positron production for observers 
close to the horizon will be dominated by in-scatter of the electrons and positrons 
into the black hole. Given coordinate dependent MDR quantities, we then 
parallel transport each of these 4-dimensional momenergy vectors parallel to 
themselves out to the inner shell of the collapsar (say at r = 50M). We then 
determine the time dilation and bin the results as a function of the poloidal 
angle. In this way we can estimate the MDR from neutrino scattering as a 
function of the angular momentum of the black hole. 

This treatment is obviously non-covariant. However, as the astrophysical 
uncertainties of the environs of such a black hole/accretion disk system appear 
to be, at present, far greater than the issues with relativity we raise here; we feel 
that our violation of general covariance to simplify the computational require- 
ments may be justified in order to provide at least a qualitative understanding 
of the off-axis energy-deposition at the inner envelope. We acknowledge that 
for a more complete characterization of the source term, a covariant treatment 
should ultimately be undertaken. In the meantime, we feel that the results 
presented herein provide a useful qualitative estimate. 

We have provided a pictorial representation of the total MDR calculated 
for a field of observers in the z-x plane perpendicular to the accretion disk in 
Fig H Here we have placed 10 4 observers in this plane and calculated the MDR 
for each using Eq. (|2^). The observers were evenly spaced in 100 equal BL 
increments (Ar) in radius from r — 2.5M out to r = 12M, and equally spaced 
in 100 BL poloidal increments (A0) from = to 9 = tt/2. The MDR densities 
were transformed to BL coordinates. The 3-volume assigned to each observer 
is, as mentioned above, coordinate dependent (AV = 2tt ^fgTj Ar A0) The last 
column of the figure illustrates the structure of the MDR generated above the 
surface of the accretion disk. 

In Fig. H we integrated the BL MDR 4-momentum out along their geodesies 
until they either hit the black hole or they reach r = 50M. We do this for 
five black hole/disk systems. The five systems correspond to specific angular 
momentums with the values, a = 0.0,0.25,0.5,0.75 and 1.0. For each value of 
a the BL time component of the MDR in the BL frame reaching r = 50M was 
binned into 20 bins ranging from 9 = to 9 = tt/2. The relative magnitude of 
these curves increase with increasing values of the specific angular momentum of 
the central black hole. We have investigated the origin of this increase and find 
that the source of the increase is the region of the accretion disk near the inner- 
most stable orbit. This is illustrated in (Fig. where we compare the angular 
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MDR profile for an accretion disk extending from 6M out to lOAf with (1) a 
Schwarzschild black hole and (2) and extreme Kerr black hole. The magnitudes 
are approximately equal and only a slight tightening of the angle of the peak 
MDR is apparent for the extreme Kerr black hole system. 

We demonstrate convergence of our results with regard to the number of 
observers we use in the z-x plane by comparison of a high and lower resolution 
run (Fig. ||). The higher resolution, which we use in our calculations in this 
manuscript, is a little smoother than the lower resolution run; however, the 
both compare well with each other. 

The relevant quantity reported by AF was the total MDR at a distant ob- 
server. If we integrate each of the five MDR curves shown in Fig. || over 17 at 
r = 50A/ we will obtain the integrated MDR into the jet as shown in Fig. ^|. It 
will be concentrated along the axis and peak along a cone of opening angle 7r/4. 
This opening angle is due to the tilt of the MDR near the surface of the disk. 
The tilt is in the direction of the disk's rotation. We find a linear dependence 
for this integrated MDR as a function of the specific angular momentum of the 
black hole. This is an extension of the earlier results because we are including 
the relatively large off-axis source term - the MDR created just above the sur- 
face of the disk via v v scattering. This disk MDR produces a 10 to 20 fold 
enhancement to the total MDR energy deposited at r — 50M over and above 
the on-axis MDR, even for a Schwarzschild black hole. Such a disk-generated 
MDR was also found in the non-gravitational models of AF.(T^| The dependence 
of this disk-driven source of momenergy on the specific angular momentum of 
the black hole is roughly linear. We predict a twofold increase in the MDR for 
an extreme Kerr black hole relative to a Schwarzschild hole. In addition, in 
(Fig. ||) we only see a slight tightening of the cone of the momenergy deposition 
with increasing a. It would be interesting to observe the jet dynamics through 
the outer envelope of the star with sources of this kind. 

7 Conclusions: Enhancement of the Energy De- 
position from a Disk-Driven Wind 

The goal of the research in the GRB central engine modeling community over 
the last few years was to look for a needed enhancement of a factor of 10 in the 
MDR. Previous work by AF and SW looked for such an enhancement over the 
non-gravitating models by examining the GR effects on neutrino pair scattering 
along the symmetry axis. They did not find this enhancement; however, in this 
paper we identify the off-axis disk-driven MDR component that yields such an 
increase over corresponding on-axis MDR values. Furthermore, this disk-driven 
momenergy occurs in the non-rotating black holes as well as the rotating black 
holes, and is even present in the non-physical models with no gravity effects (e.g. 
Fig 4 of Ref. pT|). We therefore have strong indication that the the current 
fireball central engine model discussed here is viable. Of course, more detailed 
astrophysical simulations are needed which include the complex environs of the 
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Figure 6: We bin the MDR per unit solid angle as a function of the poloidal 
angle (9) at r — 50M for black hole/accretion disk systems with five values if 
the specific angular momentum (a = 0, .25, .5, .75 and 1). Although dE/drdU 
levels off and yields results consistent with the conjecture of AF that the MDR 
is constant for small angles away from the axis; nevertheless, we do find a 
substantial disk-driven wind source off-axis. We demonstrate here the following 
two features: (1) The largest contribution to the MDR originates just above the 
disk giving factors of 10 to 20 above the on-axis non-relativistic contributions 
calculated by others. This gives rise to a disk-driven wind. This wind, when 
propagated out to r — 50M along geodesies produces the peak in the MDR 
profiles shown in the figure. This peak corresponds to depositing the bulk of 
the MDR in a cone with opening angle 9 » 7r/4; and (2) The MDR only mildly 
depends on the specific angular momentum of the disk and in an approximately 
linear fashion. The excess MDR for the extreme Kerr system results from the 
size of the accretion disk. In the Kerr system the disk extends from just above 
r = M out to r — 10M; however, in the Schwarzschild case the disk extends 
from r = 6M out to r = lOAf. These runs were performed using 10 4 observers 
located at constant Sr and 58 intervals in the z-x plane. The observers were 
placed between r = 2.5M and r = 10M. The MDR integral was calculated at 
each of these observers using 2 x 10 grid points covering the celestial sphere. 
Values of the MDR can be obtained by multiplying the numbers in the graph 
by a and by the cube of the black hole's mass. 
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Figure 7: Demonstration that the MDR is primarily driven by the size of the 
disk and thus the location of the inner-most active region (active in the sense 
of neutrino production) of the accretion disk. In the previous figure (Fig. ^) wc 
saw a factor of two increase in the peak MDR for an extreme Kerr black hole as 
compared to the Schwarzschild system. This is because the disk around the Kerr 
black hole extended down to the horizon, while in the Schwarzschild spacetime 
the disk ended at the inner most stable orbit at r = 6M. In this figure we 
compare these two systems with the same disk. In both simulations (a = and 
a = 1) we use an isothermal accretion disk with inner and outer radii of 6M 
and 10M, respectively. However, one simulation uses an extreme Kerr black 
hole (a = 1), while the second simulation uses a Schwarzschild hole (a = 0). 
The effect of angular momentum, sans location of the inner most stable orbit, 
has a marginal effect on both the magnitude of the total momenergy deposited 
at 50M and on the direction of the momenergy. Values of the MDR can be 
obtained by multiplying the numbers in the graph by a and by the cube of the 
black hole's mass. 
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Figure 8: Demonstration of the convergence of our calculations of the MDR 
at 50Af using different number of observers located in the z-x plane around an 
M = 1 Schwarzschild black hole/ accretion disk system. The high resolution 
curve uses 1.5 x 10 4 observers in the z — x quadrant (equally spaced between 
r = 2.5M and r — 12M), while the low resolution curve uses 1 x 10 4 observers 
in this quadrant. In both cases the MDR scattering integral used 100 x 200 grid 
points. Here we plot the MDR per unit 9 as opposed to MDR per unit f2 in the 
other figures. Values of the MDR can be obtained by multiplying the numbers 
in the graph by a and by the cube of the black hole's mass. 
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Figure 9: We plot the total BL time component of the MDR at 50M as a 
function of five values of the specific angular momentum (a) of a solar-mass 
black hole with a 10 MeV temperature accretion disk. In other words, we 
integrate each of the five profiles in Fig. || over the solid angle ft. This shows 
that the MDR is roughly a linear function with respect to the specific angular 
momentum (a) of the black hole. We provide a least-squares fit of the five 
numerical points. The integration of the MDR profiles are done over the entire 
northern hemisphere. If we integrated our profiles only over a small solid angle 
around the symmetry axis we would find our result are consistent with those 
reported by AF. However, we include here the disk-driven wind created by v— 
v scattering above the disk's surface. We see a large factor (approximately 
10-20 times) in the MDR from the disk surface over and above the on-axis 
source calculated by AF and SW. This enhancement occurs equally well for 
Schwarzschild black holes as for Kerr black holes. 
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black hole disk system, geometry and transient behavior of the disk and black 
hole, and a simulation of the induced plasma transport of the scattering out to 
the fireball region. Since this disk-driven MDR source is relatively large, and 
since the baryon pollution problem p9[ remains an open issue in this field, [jjj 
we argue that it may be premature to discount this off-axis source. We should 
reserve such a decision until more detailed GR, hydrodynamic, MHD and plasma 
transport calculations have been done. 

As we outlined in the paper, we have extended the earlier work of AF to 
calculate the off-axis contribution to the MDR. When using the assumptions 
made by these authors, our results are consistent their work; however, we find 
the generation of a substantial MDR above the disk, which, when added to 
the on-axis calculations of these authors produce a substantial enhancement to 
the deposition of energy-momentum. In particular, we find an approximately 
linear dependence of the total MDR at large r as a function of specific angular 
momentum of the black hole. 



dE 
d7 



soAf ylQMeVJ \M Q 



4.06 x 10 51 (1 + -73a) ( Jg*±) "(-?£- V ^ (32) 



Here we did a least-squares fit to the numerical results in Fig. |9j, multiplied by 
a times the cube of the solar mass (Mq = 1.47 x 10 5 cm). The bulk of the MDR 
at large r is deposited in a conic region centered on the symmetry axis of the 
system with an opening angle of approximately 7r/4. 

Although we provide a rigorous general relativistic calculation of the MDR 
at an arbitrary observer, computational constraints required us to make some 
approximations regarding the transport of the scattered energy out to the inner 
envelope (r = 50M) of the star. Assuming (1) the electrons produced have 
substantially higher energy than their rest mass, (2) that the electrons propagate 
away in a collisionlcss fashion and (3) that the observationally-relevant electron 
creation distribution functions are peaked around the total MDR calculated, 
then our results should be rigorously correct for a Schwarzschild black hole, and 
qualitatively correct for spinning black holes. We would not be surprised if a 
general-relativistically rigorous calculations would yield different quantitative 
results for the extreme Kerr black hole close to the horizon. However, one of 
the main results of this paper is the identification of the off-axis disk-driven 
MDR as providing an enhancement to the energy deposition rate rather than 
the GR-driven enhancements. 

We outlined in the last section an approach to calculate the electron distribu- 
tion function at an observer O at large r in a relativistically covariant way. This 
involved integrating the electron-creation distribution functions along geodesies 
into the past which originated at O. One step in this direction would be to 
extend the work here and calculate the electron creation distributions in the 
frame of our observers as opposed to the total MDR. This is not difficult in 
principle but is computationally demanding as it would require evaluating a 
6-dimensional scattering integral as opposed to the 4-dimensional integral we 
evaluated numerically in Sec. ||, Eqn. 28. The utility of such a calculation would 
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be to evaluate the point along the axis on the border between out-scatter and 
in-scatter (i.e. when the total MDR was zero). Would the electron creation 
function at and about this event yield a non-negligible contribution to the elec- 
tron distribution function at large r? 

On a more speculative note, we also are interested in calculating secondary 
scattering of the electrons and positrons within the ergosphere of the black 
hole. We have found evidence of significant in-scatter into this region. Will the 
Penrose process produce some observable high energy events from this region 
that can give us more clues to the inner workings of the central engine of these 
low mass black hole/accretion disk systems? 

One question seems of paramount importance in this field: What observa- 
tional features of a GRB can give clues as to the inner-workings of the central 
engine ? [pOf Surely future X-ray, 7-ray and gravity wave observatories comple- 
mented with theory can soon unveil these cosmic explosions. 
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